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چگیده در شرایطی که حملات تروریستی به اماکن و تأسیسات غیرنظامی روندی افزایشی پیدا کرده‌است, توجه به طراحی ساختمان‌ها در 
مقابل بارهای ضریبه‌ای اشی از انفجا از اهمیت بیشتری برحوردار است و توجهات و ملاحظات ویژه‌ای را می‌طلبد. شناسایی ماهیت بارهای 
انفجاری و نحوه‌ی اعمال آن بر سازه یکی از مسائل مهمی است که باید مورد مطالعه و توجه قرارگیرد. انفجا رآزاد شدن ناگهانی و سریع 
حجم بسیار زیادی از انرژی است که تولید نور, گرماء صدا و موج ضربه‌ای می‌نماید. این موج شامل هوای بسیار متراکمی است که به‌صورت 
شعاعی و کروی با سرعتی در حدود سرعت صوت از منبع انفجار به‌سمت خارج حرکت می‌کند. اعمال این بار هولناک بر سازه باعث ایجاد 
نغییرشکل‌های زیادی در محدوده‌ی غیرعطی و خرابی قسمت‌هایی از سازه می‌گردد که انجام تحلیل‌های غیرحطی را ضروری می‌سازد. د 
این ماله بارهای انفجاری و پارامترهای مهم انفجار معرفی شده‌اند و روابط و مدل‌های مختلف با رگذاری اشجاری مورد مطالعه و مفایسه قرار 
هل ی ماست نی ملل بر تایه بر اقهال این اوه بو بازماش بامانی ای و مرو ام اده قرار رای تب ی یا 
بهکا رگیری نرم‌افزارهای تحلیل سازه, پاسخ غیرنعطی قاب‌های مختلف ساحتمانی در مقابل بارهای ضربه‌ای ناشی از انفجار محاسبه و مورد 
ارزیابی قرا رگرفته‌است. هم‌چنین با مقایسه‌ی رفتار غیرخطی قاب‌های ساختمانی با تعداد دهانه‌ها و طبقات مختاف, عوامل تثی رگذار بر پاسخ 
سازه شناسایی و معرف یگردیده‌اند. 


واژه‌های کلیدی تحلیل دینامیکی سازه موج انفجار سازه‌های ساختمانی بارگذاری انفجاری, بارهای ضربه‌ای. 
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* تاریخ دریافت مقاله ۹۰/۲/۱۱ و تاریخ پذیرش آن ۹۲/۲/۱۰ می‌باشد. 
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(۲) نویسنده‌ی مسوژول: دانشیان دانشکده‌ی مهندسی عمران و محیط زیست. دانشگاه صنعتی شیراز. 
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معد مه 

حادثه تخریب ساختمان دولتی فدرال امریکا در ایالت 
اوکلاهاما در سال ۱۹۹۵ که ناشی از انفجار ۱۸۱۶ 
کیلوگرم 1 در فاصله ۶/۷۵ متری بود [1], موجب 
شد که بررسی انفجار و آثار مخرب آن بر ساختمان, 
مورد توجه ویژه قرار گیرد و منتج به چاپ مقالات و 
آیین‌نامه‌هایی در زمینه‌ی طراحی ساختمان‌ها و فضاهای 
شهری در مقابل انفجار گردید [2,3]. ولی آنچه در این 
آیین‌نامه‌ها به‌وضوح دیده می‌شود نوپایی و درحقیقت 
ضعف این آیین‌نامه‌ها در پاره‌ای از موارد. از جمله 
نحوه‌ی اعمال این بار بر سازه و موارد دیگری از این 
دست است. 

در ضربه‌ی ناشی از انفجار زمان بارگذاری بسیار 
کوتاه است. همان‌طور که در شکل (۱) مشاهده 
می‌شود. در طول چند صدهزارم ثانیه» فشار از فشار 
اتمسفر(و۴) به فشار بیشینه (۳:0) می‌رسد و سپس در 
اه ام ات درشاره ان تنم 
بازمی‌گردد. و در ادامه در طول چند دهم ثانیه فشار 
ناه تیش و فا ای باس | 
انفجار (۲:0)» با دور شدن از مرکز انفجار به نحو 
چشم‌گیری کاهش می‌یابد. ولی برخلاف آن» دوره‌ی 
بارگذاری (زمان اثر بار ناشی از انفجار بر سازه) با دور 


شدن از مرکز انفجار افزایش می‌یابد [4] 
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شکل ۱ تصوير شماتیک تاریخچه زمانی فشار ناشی از انفجار [4] 
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در بررسی‌های مربوط به انفجار. گام اول محاسبه 
و پیش‌بینی بار ناشی از انفجار بر روی ساختمان 
می‌باشد. برای نیل به این هدف آزمایش‌های میدانی 
زیادی صورت گرفته‌است که نتایج آن‌ها در مراجع 
متعددی آورده شده است [4,5]. هم‌چنین در اين زمینه 
نرم‌افزارهایی نظیر نرم‌افزار 000۷۷۵0 که به‌وسیله‌ی 
اداره‌ی مهندسی ارتش ایالات متحده توسعه یافته است 
[4] و یا نرم‌افزار 31854-.۸.۲ که توسط شرکت ۸۸ 
تهیه شده است [6]» موجود می‌باشد. برای محاسبه‌ی 
بار ناشی از انفجار روش‌های دیگری نیز ارائه 
شده‌است که اکثر آن‌ها بر فرمول‌های تجربی استوار 
هر اگم با تمه شا سار یط با 
(006) درسال ۱۹۵۵ ارائه گردید [7]. پس از آن 
روابط دیگری توسط سایر محققان نظیر نیومارک و 
هنسن [8] و هم‌چنین میلز [9] ارائه گردید . برای 
به‌دست آوردن زمان تداوم بارگذاری انفجار نیز روابط 
مختلفی توسط محققان ارائه گردیده است که از آن 
جمله می‌توان به روابط ارائه شده توسط مایز و اسمیت 
در سال ۱1۹۹۵ لم و همکاران در سال۲۰۰۶ و ایزدی 
فرد و ماهری در سال ۲۰۱۰ اشاره تمود [10,11,12] 

هم‌چنین در سال‌های اخیر محقفانی نظیر لم و 
همکاران تلاش‌هایی برای دست‌یابی به روش‌های 
طیفی و تهیه‌ی طیف طراحی در مورد بارگذاری 
انفجاری انجام داده‌اند [11]. در زمینه‌ی اثر انفجار بر 
سازه‌های مختلف و رفتار ساختمان‌ها در مقابل این 
بارگذاری مطالعات زیادی توسط محققان صورت 
گرفته است [13,14]. نتایج این بررسی‌ها بیانگر ایین 
مسأله است که نتایج حاصل از روش‌های تحلیلی 
اجزای محدود و نیز روش های تحلیلی که از روابط 
عددی استفاده می کنند. هم‌خوانی خوبی با نتایج 
آزمایشگاهی و میدانی دارند. در سال ۲۰۰۹ لی و 


همکاران اثر بار انفجار را بر ستون‌های فولادی با 
مقاطع مختلف بررسی نمودند [15]. در سال ۲۰۱۲ 
هادیان فرد و فراهانی تعداد مقاطع بیشتری از ستون‌ها 
را مورد ارزیابی قرار دادند و رفتار غیرخطی آن‌ها را 
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مطالعه نمودند [16]. با توجه به استفاده‌ی روزافزون از 
نرم‌افزارها برای تحلیل ساختمان‌ها در مقابل انفجان 
لازم است روش‌های مناسب برای برآورد بارهای ناشی 
از انفجار و نحوه‌ی اعمال آن بر سازه شناسایی گردد و 
دقت هر روش مورد ارزیابی قرار گیرد. این امر یکی از 
اهداف اصلی مقاله حاضر می‌باشد. 


محاسبه‌ی مشخصات مربوط به بارهای 
ناشی از انفجار 

برای محاسبه‌ی بار ناشی از انفجار روش‌های زیادی 
ارائه شده‌است. در واقع همه‌ی پارامترهای مربوط به 
انفجار وابسته به دو پارامتر مستقل هستند: یکی میزان 
مواد منفجره که با معادل‌سازی برحسب وزن 1 
معادل بیان می‌شود (1۷) و دیگری فاصله‌ی بین مرکز 
انفجار تا محل اثر موج انفجار (18). حاصل این دو 
کمیت. به‌صورت پارامتر 2 (فاصله مقیاس شده) به 
شرح زير بیان می‌شود: 


7 
در رابطه‌ی بالا 15 بر حسب متر و ۷۷ بر حسب 
گذرکوم بکار مس ژوند 119 
رابطه‌ی پیشنهادی براد (8006) برای برآورد 
فشار بیشینه‌ی ناشی از انفجار به‌صورت زیر می‌باشد 
[7]: 
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هم‌چنین رابطه‌های پیشنهادی توسط نیومارک [8] 
و میلز [9] به‌ترتیب در روابط (۳) و (4) نشان داده 


شده است. 
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(( 
اگر موج انفجار در مسیر خود به‌مانم برخورد 
کند. پدیده‌ی انعکاس موج, پیش‌خواهد آمد که میزان 
این انعکاس. متناسب با زاویه ی برخورد موجء متفاوت 
خواهد بود. در صورتی که زاویه ی برخورد صفر باشد 
(مانع موجود در مسیر» بر مسیر حرکت موج عمود 
باشد). میزان انعک‌اس ایجاد شده از رابطه‌ی زیر 


محاسبه می‌شود [13]: 
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به‌صورت رابطه‌ی لگاریتمی زیر پیشنهاد گردیده ندمت 
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هم‌چنین رابطه‌ی جدیدتر ارانه شده توسط 


ایزدی‌فرد و ماهری مطابق زير می‌باشد [12]: 
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تاریخچه زمانی فشار ناشی از انفجار که در شکل 
۱ 
فرایدلندر (۳۲۱60120067) مطابق رابطه‌ی (۸) بیان 
می‌شود [17] 
(۰ (ع ط-) دوه (عا -۲,)1 - ).۲ 
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که در آن ما زمان رسیدن موج انفجار به سازه وا 
ایکا و ۵ فرش کاهی کشت اقا کر 
ضربه ب‌صورت سطح زیر نمودار فشار- زمان می‌باشد 
(سطح زير منحنی شکل (۱) در مدت زمان یا) که 
رابطه‌ی انتگرالی آن به قرار زیر است: 
)0۱ ,۳ 1 2 .1 

در بسیاری از موارد» برای سهولت. نمودار (۲:60 
به‌صورت مثلشی فرض می‌شود (به عبارت دیگر 
تغییرات فشار با زمان خطی فرض می‌شود). که در آن 
فشار اولیه‌ی‌م:] در مدت زمان و۲ به صفر می‌رسد. بار 
انفجار به‌صورت باری گسترده بر نمای خارجی 
ساختمان اعمال می‌شود که مقدار آن با دور شدن از 
منبع انفجار کاهش می‌یابد. معمولا توزیع این بار 
به‌صورت گسترده فرض می‌شود که در ارتفاع ساختمان 
به‌صورت خطی تغییر می‌کند. البته در برخی از 
تحفیقات این بار ب‌صورت بارهای متمرکز در تراز هر 
طبقه قرار داده شده‌است. در شکل (۲) توزیع حطی و 
متمرکز بار در ارتفاع ساختمان نمایش داده شده است. 
در این مقاله. ابتدا برای انفجارهای مختلف و در 
فواصل متفاوت. بیشینه فشار و بیشینه فشار انعکاسی 
به‌وسیله‌ی نرم‌افزار ۸۰1-81256 و روابط بالا محاسبه 
شده و با هم مقایسه گردیده‌اند. پس از آن قاب‌های 
متفاوتی. هرکدام یک‌بار با بارهای متمرکز» و یک‌بار با 
بارهای گسترده‌ی خطی, با توزیم خطی در طول زمان, 
و یکبار با توزیع غیرخطی در طول زمان. تحلیل 


شده‌اند و نتایج مقایسه گردیده‌اند. 


شکل ۲ انواع بارگذاری به‌صورت خطی و متمرکز 
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وهی ال سای ی تخل سازهای شمان 
مختلف در نرم‌افزار 
بارهای انفجاری به‌دلیل داشتن طبیعتی ضربه‌ای» در 
زمان بسیار کوتاهی به سازه اعمال می‌شوند و به همین 
دیل با پارهای جانبی دیگر تظیر باد و زلرله متفاوت 
کی از الا اس رس بانهای فامس از 
انفجار ابتدا برای حدود پنج هزار موقعیت متفاوت در 
فواصل ۲ تا ۸۵ متری و حجم مواد منفجره از ۱۰ تا 
۰ کیلوگرمی, فرمول‌های تجربی ذکر شده در بخش 
قبلی مقاله. محاسبه شده و نتایج حاصل از آن‌ها با 
نتایج حاصله از نرم‌افزار 31850-.۸.1 مقایسه 
گردیله‌است. 
در مرحله‌ی بعد. سازه‌هایی که به‌صورت کامل و 
با رعایت تمامی ضوابط آیین‌نامه‌ای و برای بارهای 
لرزه‌ای طراحی گردیده‌اند. مورد ارزیابی قرار گرفته‌اند. 
بدین منظور قاب‌های سازه به‌صورت مجزا تحت بار 
ناشی از انفجار قرار داده شده‌اند و با بررسی نتایج 
نکاتی درباره این نوع بارگذاری بیان گردیده‌است. 
سازه‌های ساختمانی مورد مطالعه در شش حالت 
متفاوت. در دو. پنج و هشت طبقه و هر کدام یکبار با 
تاه از کیان اش تاه تسا اه امه 
ساختمان‌ها با اهمیت زیاد و در منطقه‌ای با خطر 
زلزله‌ی بسیار زیاد. قرار داشته وبا استفاده از فولاد نرمه 
ماتفماین با شش ری فان ۱۱۳9 295 عت زیت 
پواسون ۱0.3 و مدول الاستیسیته 0۳۵ ۳-200 
به‌صورت قاب خمشی ویژه‌ی فولادی طراحی 
دنله ال تفر شک ۳ وپمات قاب‌ها نشاق 
داده شده‌است. 
فرضیات به‌کار گرفته شده در تحلیل‌ها عبارتند از: 
هاش ات اقفر یه جالی سارت مها وه 
اعمال شده‌است؛ در حالت اول, بارگذاری با توزیع 
گسترده بر روی وجه ساختمان ( وجه ستون) و از 
نظر تاریخچه زمانی مطابق با شرایط واقعی در نظر 
گرفته شده‌است. در حالت دیگر بار انفجاری 
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به‌صورت گسترده بر روی وجه ساختمان ولی با 
توزیع مثلشی در طول زمان و در حالت سوم 
به‌صورت بار متمرکز از نظر تاریخچه زمانی مطابق 
با شرایط واقعی در تراز طبقات در نظر گرفته 
شده‌است. با توجه به جدول موجود در 
8 مراد منفجره یک خودرو سواری در 
حدود ۳۵ تا ۸۰ کیلوگرم خواهد بود که اين بار در 
فاصله‌ای در حدود ۱۵ متسر باعت تخریب 
میتوزن‌ه ای سا شمان گر اه شنت 181 در 
مدل‌سازی‌های صورت گرفته‌ی اولیه این نتیجه 
تأیید شد و برای مواد منفجره در فواصل کمتر از 
۰ متر. ستون‌های بیرونی کاملاً تخریب گردیدند. 
با توجه به این نکته و برای این که بتوان اثر نحوه‌ی 
اعمال بار ناشی از انفجار بر سازه را بررسی نمود و 
نیز به‌منظور بررسی رفتار غیرخطی سازه قبل از 
خرابی. در اين مقاله بار مبنای طراحی. بار ناشی از 
انفجار یک خودرو با ۵۵ کیلوگرم 13 در فاصله 
۰ متری در نظر گرفته‌شده است. نتایج حاصل از 
این تحلیل برای سازه‌های تحت اثر انفجارهای با 
ای سس ار اسان ( در دود دا و 


متر) دارای کاربرد می‌باشد. 


۲- بار موثر بر هر قاب. با توجه به عرض باربر جانبی 


و 


عِ- 


۵ 


۱۳۵ 


آن محاسبه شده‌است که عبارت است از مجموع 
نصف فاصله آن قاب تا قاب‌های کناری. هم‌چنین 
فرض شده‌است تمام بار ناشی از انفجار از طریق 
دیوارها به سازه منتقل گردیله‌است.: 

نرم‌افزار مورد استفاده برای تحلیل قاب‌ها؛ نرم‌افزار 
5 می‌باشد و تحلیل‌های صورت پذیرفته در 
آن به‌صورت دینامیکی تاریخچه زمانی و با درنظر 
گرفتن رفتار ارتجاعی- خمیری کامل برای فولاد و 
با استفاده از جعبه‌ی ۳ این نرم‌افزار انجام 
گرفته‌است. 

المان مورد استفاده در تحلیل 3۸1۷189 می‌باشد 
که المانی با امکان رفتار غیرخحطی. درجه دو و با 
شش درجه آزادی در هر گره به صورت سه درجه 
آزادی جابه‌جایی و سه درجه آزادی دورانی است 
[19]. 

مقاطع مورد استفاده در ستون‌ها برای قاب دوطبقه: 
0 1۳ و طبقات یک 0 ۳و برای قاب 
هشت طبقه: طبقات اول تا سوم 0 ۰.1۳ طبقات 
چهارم تا ششم 360 1۳3 و دو طبقه‌ی آحر 1۳8 
0 است. هم‌چنین مقاطع تیرها شامل 1۳1270 و 
0 می‌باشند. 


له و له و مه 4 صلم4 ماه صله 4 صاه 4 صاه 4 ما 


سال بیست و ششم. شماره یک. ۱۳۹۳ 


نشریه مهندسی عمرآن فردوسی 


۱۳۹ 


نتایج حاصل از روش‌های محاسبه‌ی بار انفجار و 
مقایسه‌ی آن‌ها 

با توجه به این‌که نتایج نرم‌افزار ۸۰۲-3185 در 
حقیقت حاصل از داده‌های میدانی و آزمایشگاهی 
می‌باشد و مورد تأیید مراجع معتبر علمی و نظامی 
است. در بسیاری از مقالات به‌عنوان معیاری برای 
کنترل و بررسی نتایج محاسبات مربوط به میزان بار 
انفجاری مورد استفاده قرار گرفته‌است [15]. به همین 
دلیل در این بررسی نیز از نتایج حاصل از اين نرم‌افزار 
برای کنترل دقت روابط تجربی استفاده شده است. در 
ابتدا بارهای انفجاری حاصل از مواد منفجره با 
وزن‌های مابین ۱۰ تا ۲۰۰۰ کیلوگرم 1۲ (۸۰ ۲۰ 
۳۰ ۰۰ کیلوگرم) و فواصل مختلفی بین ۲ تا 
۵ متر (۲ ۳ ۶. ...تا ۸۵ متر)» توسط نرم‌افزار 
۸.1.2 به‌دست آمده و سپس نتایج حاصل از 
فرمول‌های تجربی نیز محاسبه گردیده‌اند و اين نتایج با 
هم مقایسه شده‌اند. 

در مجموعه شکل‌های شماره‌ی (۶) نتایج حاصل 
از محاسبه‌ی زمان تداوم بارگذاری ,4 برای وزن‌های 
۵ ۵۱۰18۵ 12 ۱۱۰۰ و 12 ۲۰۰۰ مواد منفجره 
1 که از فرمول‌های تجربی پیشنهادی توسط 
ایزدی‌فرد و ماهری (۱۷ »8 1) و تم (1870) محاسبه 
شده‌اند. با نتایج حاصل از نرم‌افزار ۸۰1-8189 مقایسه 
و میزان حطا نمایش داده شده است. 

چنانچه از این شکل‌ها پیداست. رابطه‌ی مربوط 
به ایزدی‌فرد و ماهری در هیچ فاصله و با هیچ میزان 
خرج انفجاری. خطایی بیش از ۵ درصد ایجاد نمی کند. 
یرس رید وس یش 1 
حدود ۳۰ متر و برای وزن‌های مواد منفجره کمتر از 
۹ کیلوگرم 1 درصد خطای این رابطه بسیار 
کم خواهد بود تا جایی که در فواصل بیشتر به صفر 
نزدیک می‌گردد. ولی درباره‌ی رابطه‌ی مربوط به لم 
باید گفت که اين رابطه در فواصلی کمتر از حدود۱۰ 
متر دارای دقت مناسبی نمی‌باشد و درصد خطای آن تا 
حدود ۱۵ درصد هم می‌رسد که اين میزان خطا نسبتا 


زیاد است. ولی این رابطه در حجم‌های مواد منفجره‌ی 
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و ان ۳۵۳ هبل رتیت در فراع شیق ۱۸ 
تا ۳۰ متر از دقت بهتری نسبت به رابطه‌ی ایزدی‌فرد و 
ماهری برخوردار است ولی در فواصلی غیر از این 
میزان خطای این رابطه به‌مراتب نسبت به خطای 
رابطه‌ی ایزدی‌فرد و ماهری بیشتر خواهد بود. همین 
وزن‌های معمول مواد منفجره. هم‌چنین محدود بودن 
میزان خطای رابطه‌ی ایزدی‌فرد و ماهری به کمتر از ۵ 
درصد. باعث ایجاد اطمینان بیشتری به این رابطه نسبت 
به رابطه‌ی لم می‌شود. 

در مجموعه شکل‌های شماره‌ی (۵)؛ نتایجم حاصل 
از مقایسه‌ی بیشینه فشار واقع ۳۹0 برای حجم‌های 
8۵ ۶۰ 12 ۵۱۰ 12 ۱۱۰۰ و 12 ۲۰۰۰ مواد منفجره 
۲ بین نتایج حاصل از فرمول‌های تجربی پیشنهادی 
مشاهده می‌شود. 
به براد در غالب فاصله‌ها و با غالب میزان خرج‌های 
هم‌چنین کاملاً مشاهده می‌گردد که در فواصل بیش از 
حدود ۸ متر درصد خطای این رابطه بسیار کم (کم‌تر 
از ۰/۵ درصد) خواهد بود و در فواصل بیشتر به صفر 
نزدیک می‌گردد. ولی درباره‌ی رابطه‌ی مربوط به میلز 
باید گفت که این رابطه در فواصل کمتر از حدود ۸ متر 
دارای درصد خطایی حتی تا حدود ۱٩‏ درصد می‌باشد 
که این میزان خطا نسبتا زیاد است و لذا در اینن 
محدوده. این رابطه دارای دقت مناسبی نمی‌باشد. ولی 
همین رابطه در غالب وزن‌های مواد منفجره در فواصل 
بیش از ۱۰ متر از دفت به‌مراتب بهتری نسبت به 
درصد). در مجموع و با دقت به همه‌ی این شرایط 
یعنی حطای کمتر در مورد هر میزان فاصله (با تأکید بر 
فاصله‌ی بیش از ۱۰ متر) و هر حجم مواد منفجره 
استفاده از رابطه‌ی میلز نسبت به رابطه‌ی براد از دقت 
استفاده از رابطه‌ی براد نتایج دقیق‌تری به همراه خواهد 


تیال بیست و ششم» شماره یک. ۱۳۹۳ 


داشت. منفجره بین نتایج حاصل از فرمول‌های تجربی 
در مجموعه شک های شماره‌ی (1) » نتایج پیشنهادی براد و میلز با نتایج حاصل از نرم‌افزار -.۸ 
حاصل از مقایسه‌ی بیشینه فشار انعکاسی ۳۰ برای 1251 مشاهده می‌شود. 
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شکل 1 درصد خطای موجود در مقایسه‌ی روش‌های براد و میلز 
با ا8125-.۸2.3 برای محاسبه‌ی بیشینه فشار انعکاسی ۳۰ 


رابطه‌ی مربوط به براد در بیشتر فاصله‌ها و با 
میزان خرج‌های انفجاری مختلف. کم‌تر از ۶ درصد 
تحظا اناد میک هی فواصان یش زوسن 
درصد خطای این رابطه بسیار کم (کم‌تر از یک درصد) 


خواهد بود. درحالی که درباره‌ی رابطه‌ی مربوط به میلز 
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باید گفت که این رابطه دارای خطای زیادی در فواصل 
کمتر از حدود ۸متر است که این خطا تا حدود ۱۳ 
درصد می‌باشد. هم‌چنین این رابطه در همه‌ی وزن‌های 
مواد منفجره. در همه‌ی فواصل از دقت به‌مراتب 
کمتری نسبت به رابطه‌ی براد برخوردار است. در 
مجموع. برخلاف نتایج حاصل از مقایسه‌ی بیشینه فشار 
موم برای محاسبه‌ی بیشینه فشار انعکاسی :۳ با توجه 
به خطای کمتر استفاده از رابطه‌ی براد نسبت به 


رابطه‌ی میلز از دقت بهتری برخوردار است. 


محاسبه و مقایسه‌ی پاسخ قاب‌ها به پار انفجاری 
در بخش قبل روابط موجود برای محاسبه‌ی میزان بار و 
زمان اعمال آن بر سازه مورد بررسی قرار گرفتند. در 
این بخش به بررسی نحوه‌ی اعمال این بار بر سازه 
می‌پردازيم. نکته‌ی قابل توجه این است که با 
وجودی‌که دیوارهای ساختمان پس از برخورد موج 
انفجار تخریب می‌شوند. ولی با توجه به دوره‌ی بسیار 
کوتاه بارگذاری این بار عملاً قبل از تخریب دیوارها 
به عناصر باریر سازه منتقل خواهد شد. لذا در این مقاله 
با وجود این‌که دیوارها مدل نشده‌اند. ولی بار انفجاری 
زارد بر آن‌هسابه منتون‌ها فتقل شنده است.برای 
محاسبه و اعمال بار انفجاری موثر بر ساختمان از نتایج 
پیشینه بار انعکاسی مربوط به نرم‌افزار ۸.۲.8126 
استفاده شده‌است. نمودار بارهای وارد بر سازه در این 
مقاله» در شکل (۷) مشاهده می‌گردد. برای این که بتوان 
توزیع بار در زمان را ب‌صورت خطی فرض کرد این 
نمودار باید به‌گونه‌ای در نظر گرفته‌شود که فاکتور 
ضربه برای این حالت با نمودار اصلی مساوی باشد. 
چنان که در نمودار بارگذاری (شکل شماره ۷) مشاهده 
می‌شود. زمان رسیدن موج انفجاری به سازه ۰/۰۷۳۹۷ 
ثانیه است که در هنگام تحلیل سازه‌ها؛ این زمان به 
عنوان مبداً زمانی» (زمان صفر) فرض گردیده‌است. 

در شکل شماره‌ی (۸) و جدول شماره‌ی (۱)؛ 
بیشینه جابه‌جایی مرکز جرم بام در حالت‌های مختلف 
باررگذاری با هم مقایسه گردیده‌است. نتایج نشان 


سال بیست و ششم. شماره یک ۱۳۹۳ 


وس کی تن انا تقو یه یی کته کتاماه نتایج ناشی از بارگذاری گسترده. مشابهت زیادی با 
محافظه کارانه است. به‌جای بارگذاری انفجاری با نتایج حاصل از بارگذاری به‌صورت متمرکز در تراز 
توزیع غیرخطی نسبت به زمان, از بارگذاری انفجاری طبقات دارد. اما این مطابقت با افزایش تعداد طبقات 
با توزیع خطی نسبت به زمان استفاده نمود. هم‌چنین کمتر می‌شود. 
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نمودار با ر گذاری نسبت به زمان 


با رگذاری با توزیع غیر خطی در زمان 0 
با رگذاری با توزیع خطی در زمان -س- سس لك 
0 طّ 

40 

-10 

89 0.0764 0.0759 .۰ 00754 ۰ 0.0749 ۰ 0.0744 ۰ 0.0739 
زمان (ثانیه) 
شکل ۷ نمودار بارگذاری‌های مورد استفاده در اين مقاله 
بار گسترده یکنواخت غیرخطی سح _بارمتمرکز < سس _بار گسترده یکنواخت خطی 028 


شکل ۸ مقایسه‌ی جابه‌جایی تراز بام در ساختمان ۵ طبقه با توجه به شرایط بارگذاری‌های متفاوت 


جدول ۱ بیشینه جابه‌جایی مرکز جرم بام در قاب‌هاء تحت اثر بارهای متفاوت (60) 


بار گسترده با توزیع غیر خطی | بار گسترده با توزیع خطی | بار متمرکز, با توزیع حطی 
در زمان در زمان در زمان 
قاب ۲ طبقه ۸/۰ ۸/4 ۸/۹۹ 
قاب ۵ طبقه ۳۳/۷ 2 ۱۳/۱ 
قاب ۸ طبقه ۳/۱ ۳۱/۰ ۳۳/۹۶ 


نیروی برشی در ستون اول - سب سس 


نیروی برشی در ستون دوم 


زمان (دهم انیه) 


شکل ٩‏ مقایسه‌ی نیروی برشی در ستون‌های ردیف اول و دوم قاب ۵ طبقّه 


سال بیست و ششم شماره یک ۱۳۹۳ نشریه مهندسی عمران فردوسی 


با توجه به نتایج حاصل از تحلیل‌های صورت 
گرفته. به‌طور قطع می‌توان گفت که بیشترین تنش در 
ستون‌های اول سازه (نزدیک‌ترین ستون‌هابه محل 
انفجار) 23 می‌دهد. این تنش‌ها شامل تنش‌های فشاری. 
کششی و برشی هستند. هم‌چنان که در شکل‌های 
ای ی( تین اس اک توت 
بزرگ ترین نیروی وارد بر ستون‌هاء نیروی برشی است. 
تحلیل دینامیکی تاریخچه زمانی برای بارهای انفجار 
بسیار مشابه با بارهای زلزله می‌باشد. فقط با توجه به کم 
بودن زمان تداوم بار لازم است که طول گام‌های زمانی 
بسیار کوچک انتخاب گردد و در ضمن با توجه به 
آهنگ سریع کرنش لازم است که مقاومت مواد در 
ضریب افزاینده ضرب گردد [20]. 

در جدول شماره‌ی (۲) حداکثر نیروی برشی در 
ستون‌های طبقهی اول قاب‌های مختلف ارائه 
گزذیله‌است. نکعه‌ی عهم این است کنه مییران خذاکتر 
نیروی برشی ستون اول (نزدیک‌ترین ستون به محبل 
انفجار) با افزايش تعداد دهانه‌ها تغییر چندانی نمی‌یابد. 


نیروی محوری در ستون سوم 


نیروی محوری در ستون دوم 


بارگذاری انفجار و ارزیابی پاسخ غیرعطی سازه‌های .. 


با توجه به این‌که اين بیشینه نیروی برشی در همان 
زمان‌های اولیه (حدود ۰/۰۰۰۵ تا ۰/۰۰۵ ثانیه آاغازین 
بارگذاری) رخ می‌دهد و هنوز کلیت سازه پاسخی به 
بارگذاری نداده است. مشاهده می‌شود که این میزان 
پیشینه در ستون طبقه‌ی اول با افزايش تعداد طبقات 
ساختمان افزایش چشم‌گیری نمی‌یابد. هم‌چنین مشاهده 
می‌شود که بیشینه نیروی برشی در دیگر ستون‌ها باهم 
برابر است. 

در شکل شماره‌ی (۱۱) جابه‌جایی مرکز جرم بام 
برای قاب ۸ طبقه در دو حالت سه‌دهانه و شش‌دهانه 
با هم مقایسه شده‌اند. نتایج حاصل از تحلیل‌ها نشان 
می‌دهد که قاب‌هایی که تعداد دهانه‌های بیشتری دارند. 
جابه‌جایی کلی کمتری نسبت به قاب‌های با دهانه‌های 
کر دارشد که این مسالهمی ترانت هدلیل ستصی. بیقر 
جانبی این قاب‌ها باشد. به عبارت دیگر با افزایش 
تعداد دهانه‌ها (تعداد ستون‌ها در یک قاب عمود بر 
جهت انفجار) می‌توان جابه‌جایی‌های ناشی از انفجار را 
در ساختمان کاهش داد. 


نیروی محوری در ستون اول - س - < 


۳ ۱ ۳ 30 
۱ ِ !۱" ۱ 
۱ ۱ ۱ ۱ 
متسه 1 20 
٩ ۱ ۱ ۱ ۱ ۱ ۱ ۱‏ 
ات ۱۱۱۱۱ ۲ 3 
1 [ ۱۱ ۱۱ ۱ ۱۳ ۲۱ 1 اما ۲ وا ۰ 1 ۳1 [۳ ۸ 0 « 
٩ ۱۱ ۶ ۲ ۶ ۱ ۱۲۱۱۲ ۱ ۱ ۱ 1 ۱۹ ۰ ۱‏ 1 
ب ۱ ,1 / ۱ با ما ۷ ۱۱ ا اما اب ام ما اه ۷ 40 
جا1 ۶3۸ ۵ ما۱ ۱ ۹۹ ۳ 4 ۵ 2 ۹41 مرجم( ام( ۸ ۵ ۳-۳۸۵۳ 9 
مت 
۱ ل(( ۷ ۷ 1 ۱ ۱۱ ۱۱۱ ۷ ۲۱۲ ۲ را / 909 3 
4 1 ۱ و ۳ ۷ با 1 ۷ 3 
1 ۱ , ۱۱ ۱۱ ۱۱ ۱ ار با ۲ ۱ ۳ ۱ ۸ ۱۳ : 
سروس اس سس سس تیا + ۷ 0 « 
۸ نا ۱ از نا ۱ ۱ ۱ را زا را ۱۱ ۲ ۷۱ تن 
تا ۱ ۱ ۱۱۱ ۷ 260 
۱ 3 15 ۳ ۲ 
زمان (دهم انیه) 
0- 
شکل ۱۰ مقایسه‌ی نیروی محوری در ستون‌های ردیف اول» دوم و سوم از قاب ۲ طبقه ۳ دهانه 
جدول ۲ حداکثر نیروی برشی ستون‌های طبقه‌ی اول. در قاب‌های متفاوت (کیلونیوتن) 
۱ ۱ مه | ۳ | مق ىا ۱ را ۷ 
ی ۱ ۱ 
دهانه دهانه دهانه دهانه دهانه دهانه 
نی خشر ان ۱9۲۵ ۱9۲۹ ۱۸۳۸ ۱۸۳۸ ۳۳۹۳ ۳۳:۶ 
ستون ردیف دوم ۱۶ 3۳ ۱۲ ۹۸ ۳۳۲ ۱:۱ 
ستون ردیف سوم ۱ ۳ ۱۵۲ ۹۸ ۲۳۲ ۱:۱ 


نشریه مهندسی عمرآن فردوسی 


سال بیست و ششم شماره یک ۱۳۹۳ 


احمد فراهانی- محمدعلی هادیان‌فرد 


۱۳۱ 


دهانه 3قان با تست 


دهانه ۵قاب با 


032 
016 
7 
15ع6-: 
0 
6 و 
32- 


شکل ۱ مقایسه‌ی جابه‌جایی مرکز جرم بام برای قاب ۸ طبقه». در دو حالت ۳ دهانه و 1 دهانه 


برای محاسبه‌ی فشار ناشی از انفجار روابط متعددی 
ات ایا را 
می‌توان به روابط ایزدی‌فرد و ماهری و لم برای 
محاسبه‌ی دوره بارگذاری و روابط براد و میلز برای 
محاسبه‌ی بیشینه فشار اشاره نمود. با توجه به 
بررسی‌های انجام شده در اين مقاله. موارد زیر به‌دست 

امده است: 

۱- روابط ایزدی‌فرد و ماهری به‌طور کلی نسبت به 
روابط لم دقت بیشتری دارند. خطای این روابط 
حداکثر به ۵ درصد محدود می‌شود. 

۲- در فواصل کم از محل انفجار درصد خطای روابط 
م زیاد است. 

۳- در مورد روابط مربوط به بیشینه فشار ,3 در 
مجموع استفاده از رابطه‌ی میلز دقت بهتری نشان 
می‌دهند. 

6- در مورد روابط مربوط به بیشینه فشار انعکاسی :۰۳ 
بر خلاف روابط مربوط به بیشینه فشار وی استفاده 
از رابطه‌ی براده در مجموع دقت بیشتری را نشان 
می‌دهد. 

پس از مدل‌سازی سازه و بارگذاری انفجاری بر 
روی آن در نرم‌افزان نتایج زیر قابل ذکر است: 
۱- برای مدل کردن بارهای ناشی از انفجار بر روی 


روت 


سازه استفاده از بارهای متمرکز در تراز طبقات با 
و میزان حطای آن نسبت به حالت بار گسترده با 
۲- بارهای ناشی از انفجار در ستون‌های نزدیک به 
ق بان نادزی وزی رشن فا یبای 
می‌کند که این نیرو در ستون‌های ردیف‌های بعد 
۳- افزایش تعداد دهانه‌های قاب تأثیر چندانی بر 
نیروهای ایجاد شده در ستون‌های ردیف اول 
(نزدیک به محل انقجار) ندارد, اما باعث کاهش 
تغییرمکان جانبی کلی قاب می‌گردد. از طرف دیگر 
با افزایش تعداد طبقات ساختمان, نیروهای ایجاد 
شده در ستون‌های طبقه‌ی پایین ساختمان افزایش 
لت یر تم تایه ها بان هیا کر هه 
به‌دلیل بالا بودن نیروهای ناشی از انفجار در آن‌ها 
ندسبت به نبروهای ظرانخسی معمولی سازه. در 
مقایسه با ساختمان‌های بلند دارای آسیب‌پذیری 
بیشتری در مقابل بارهای ناشی از انفجار خارج از 
ساختمان می‌باشند. البته در این مورد نیاز به 
بررسی‌های بیشتر و بر روی تعداد قاب‌های 
بیشتری می‌باشد تا بتوان به‌طور قطع اظهار نظر 


نمود. 


6 ر م۸۱ طحمظ امه ۷ )فصلجعه ففصتلانباظ فمتامعا۱۳۱۵ ویزنک ۱۷۵۱62۵۳5 4ص بخ ر«مصعصم ...1 


(1996) ,51-54 .00 و و1 ,۷۵۱ مطمتامناتعومن مج صملوه‌با آمتتمنتتو من ۳۲۳۵۲1۱001621 


عامحتاظ ]19۳018 اهتاجها۳۵ ماهع۱0۵ ما انامه عمصمععلم؟ لمع صقامصما! ۵۶ امممنجمونا ولا ...2 


2003(۰) ,426 ۲۴۱۷۸ ,۸00 اجمصم‌مممه۱ رمصمع م۴ ۲۵0۵۲21 و وعصتلاتباظ )نع 


سال بیست و ششم شماره یک ۱۳۹۳ 


نشریه مهندسی عمرآن فردوسی 


۱۳۲ بارگذاری انفجار و ارزیابی پاسخ غی رحعطی سازه‌های ... 


مهع۱0 ما وعصتلاتبظ آهله‌ممصمصمن ۵ مادهنا زمر م۳ راتسمعه مصهاممم۲ ۵ کممصنجمونز وت 
2003(۰) ,427 ۲۳۱۸ ,0صوع۸ اممصه‌عمصه]۱ زمصمع تمه اح0۵ع۳ رهاظ ۵۳۳0۲۲8[ 

۶ نما ول رد عطممهع۲۷۷ آمممتامه من ۵۲ ملوهنا ۳۲۵۲۵۵0۷۵ ۵۶ قامامم‌صهلمتا۳۳ ,۷]5-855-1(] 
1990(۰) .۲.0 طماعصتطاوه ۷۷ ,تحصنض مط 

فص 0۶ کاممصانهمه۱۳ تصامز رعممادمامدظ تماجمل۱هع۸ ۵0۴ ۳۲6۵06 معط )وروم1 ما ومستتمناا۵" ,5-1300 1۷۲[ 
1990(۰) و.).(1 طماعصاطوه ۲ ,بردا صرح ۳۵۲۵۵ بخ ررض 

2007(۰) بلاعصامن عطصامها0 ۳۲۵۲۵۵۲0۷۵ رصصهتع0۳0 1261ظ-.1 :۰۸ و.عصا ر(هگه) 0012465ووظ )مومع ۵1160ج ۸ 
رز کر ۱۱۱۶۱۰۳/۸ 

,۷۵۱26 رکعاکزاظ 60:آمم( ۵0 آ۵باول ۷۵۷6۵ اعهاها لمتمطمو ۵۶ فصمتنامی لهمتمصصتاط ونان ر006ظ 
(1955) ,766-775 00۰ ,6 :195116 

ممه10 ۷ مه امک رقمتنتمنتاد اصحاعاعی؟! )فقاظ ۵۶ صمنوهنا ر.لیظ مطهکصه۲۱ مه .۱ توا 
1961(۰) ,فقو لا ر1ع۷۲۵ ۱۱6۷۲ بأان-۱۷]۵۵۲۵۷ ,۵0 0ص نت۲۱2 .ع۲ رد ,۷۵1 ۲۱۵۵00006 

"۳۲/۵66 جممرمع/۷ مصح عصمنومامد تعلعیفک؟ ما مامتطگ ماعمممن ۵۶ طعلوعن[ مط۳1 یهت وال 
-61 00۰ یکالا مطعتااظ رکط۵ 2۳۵۵6۵ ۳۵2۵۲۵ 10۲ 60۳8۵۲6۵ 0۴8 0۳/6۲۵۸۵۵) ۳ 156 ۱6 ]0 کو0 ۳۲0066 


73, )1987(, 


4 ۱۵0008 رانا 6۵۳۷16۵ ۲6۵1۶0۳0 وقحصمط 1 و کمصالاتبدظ وم فامع]۶ تفقاظ" رربایظ مطاتصه 20 .0 رونره]۷ . 


41۳. 1210., )1995(, 


6 ۱۱۵201۳82۰ )فقاظ ۶۵۲ فممتنامی نامع فقممم‌ومک وربآ رمعلا 200 .۲ رعللمه رلط محطقا . 


(2004) ,28-44 .وم رگ ۷۵۱ ,6۵و هام50 ]0 [۲۳۵اص 


لمستمنتنه ۵۶ [۷۵ع1 مطا ووعععظ 4۵ ۱۲۵/۵۵0 اند ۵۶ ممتاهلاممه ریا رته‌تطه]۱ مضه هیک ر2201]27] . 


20068(۰) وب«روما166۳00 ۵000 5668۱66 ۷۵۵۳۵۴۵6۵1 ]0 [۲۵ا0ل ,"1۵205 فحاظ مه مبارا معجمصور1 


حبذ ومتیاامنتاگ چم فامع۲۶۲6 )فقاظ مج عصتلعما تعقاظ و.ط روحعحصف فصح یه مقاموتات ر.ظ رکل۷]۵8 و. 1 رمعل . 


2007(۰) ر۱۱66۳۱۵ 9 هام5 ]۵0 0۵ 600۳0۳۱6 0۷۵۱۵۲ 


9۳۱۵۱۱۲۵ م0077 صرح 4صج منالماع)۱۷ ۶ )عجمرحضا جرج )فحاظ ۵۶ 5۵۷1۵7 ۲۸ و.ت) رتاناً 2040 .۲ ما2 . 


2007(۰) ,5۱9۱66۲۱ آه 50 ]0 0۱0۲۳۵1 


۲۵۵0108۳ )6ه۱اظ تعهصنا فصصصتامن اما ۷۷۰۵02۵60 ۵۶ وقمممووف لممم[" ر.ل تک 240 .1 رحطتکا ریک رب . 


٩0.1 00. 25-38, )2009(‏ .31 ۷۵۱ رکمتصهط۱۷]60 4صه عصاهمصتم‌ص ناه 


آمصمتاهعز وومه فصصصنتامی آهعاو ۵۶۴ ]6۲60 فطا م۳۵ .یه تحطمتطول 0صع بظ رتصفطهنه۲ ,ما رلتهگصهلله] ۲ . 


۵ ۱9۱۱۵6۳۱۱9 5۵ ]0 ۵1اه عصتل1۱۵2 اقداها م4 ۲۵4مو میاه ممطه قمتامدتتقطهها فص مرم فقت)010۵۵۲ 


.(2012) ,449-463 .0 ,4 ,0 ,44 ,۷۵1 رکی ۷۲66۳۲۵۱[ 


(1973) ,126 مطلاکتنظ ر۳۲655 1۵625 ۵ الوزمتصرن] رد ما فصمتومای<۳ رب ۷۷ رتععلوظ . 


۶۴ 0290 م1 دامع[۳:0 ٩0001‏ مهو صفزوهنز عم تم۳۳ رفظ پمممصم‌مممها مصوع‌مص لقتم0ع۳ . 


2003(۰) ,428 ۳۳۱۷۲۸ ,"ععاعمتاظ افادمه ]1 


۰ اصمصصم ۱ رممممتمعم1 اممحصماظ قامعا( ممممعع]م۳ , :11 ۸۲5۷5 . 


۴ 09۵ "5۵/69 من مرن که ومتانه۳۵۵ظ اهمنصعطامع]۱۷ 9۵۵" ر.ظیک رم 20 .ظ مصفهتطعتام۵. 
.(1987) ,663-672 .00 و4 .0 و113 ,۷۵۱ رفظ ,۱9۱۳66۲۱۵ 50۳6۵ 


نشریه مهندسی عمران فردوسی سال بیست و ششم شماره یک ۱۳۹۳ 


